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A.  Introduction  1 /jr; 

This  report  serves  to  describe  the  research  accomplished  under 
the  support  of  ARO  Grant  DAHC04-75-G-0170.  Our  efforts  have  focused  in  five 
different  areas,  all  related  to  a better  understanding  of  solid  surfaces: 

(1)  The  analysis  of  surface  structures  by  low-energy  electron  diffraction 
(LEEO)  using  both  full  multiple  scattering  (dynamic)  and  transform-convolution 
(kinematic)  numerical  schemes;  (2)  Other  research  related  to  LEED  such  as 
an  improved  method  of  measuring  experimentally  the  angle  of  incidence  of  the 
electron  beam  in  a LEED  apparatus,  the  measurement  of  surface  and  bulk  Debye 
temperatures,  and  the  evaluation  of  adparticle-adparticle  interaction  energies 
via  a measurement  of  order-disorder  phenomena  using  LEED;  (3)  An  examination 
of  desorption  phenomena  including  a theoretical  treatment  of  thermal  desorp- 
tion mass  spectrometry  (TDS)  and  experimental  measurements  of  electron 
stimulated  desorption  (ESD);  (4)  Model  calculations  concerning  the  dependence 
of  the  spectral  intensity  upon  adsorbate  concentration  in  inelastic  electron 
tunneling  spectroscopy  (IETS),  the  change  in  electronic  free  energy  upon 
reconstruction  of  an  ideal  (lxl)  surface  structure  to  either  a reconstructed 
(2x1)  or  c(2x2)  structure,  and  the  change  in  the  electronic  density  of  states 
when  a two-level  adsorbate  interacts  with  the  s-band  of  the  (100)  surface  of 
an  fee  metal;  and  (5)  Chemisorption  and  surface  reactions,  in  particular, 

Hg,  O2  and  CO  on  the  (111)  surface  of  Rh,  and  CO  on  the  (110)  surface  of  Ir. 
Each  of  these  topics  will  be  summarized  in  more  detail  in  sections  B through 
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F below.  References  cited  in  these  sections  are  to  the  publications  which 
have  resulted  from  this  grant  support.  For  convenience,  these  references 
are  listed  separately  in  section  6. 


B.  Analysis  of  Surface  Structures  by  LEED 

The  major  thrust  of  the  research  supported  by  this  grant  during  the 
past  three  years  has  been  concerned  with  a determination  of  the  structure  of 
solid  surfaces.  This  includes  both  the  relaxation  (if  present)  of  clean 
metallic  surfaces  as  well  as  the  position  of  ordered  adatoms  on  metallic 
surfaces.  This  research  has  been  both  of  a theoretical  and  an  experimental 
nature.  On  the  theoretical  side,  there  have  been  two  major  advances,  namely, 
the  development  of  a transform-convolution  method  of  analysis  of  the  structure 
of  clean  surfaces  via  the  direct  analysis  of  LEED  intensity-voltage  (I-V) 
beam  profiles,  and  the  ability  to  perform  full  multiple  scattering  (dynamic) 
calculations  for  both  clean  surfaces  and  surfaces  on  which  there  are  ordered 
overlayers.  The  latter  ability  is  due  entirely  to  Dr.  Michel  Van  Hove  who 
spent  one  year  in  our  research  group  at  Caltech  as  a postdoctoral  fellow 
prior  to  moving  to  the  University  of  Munich.  On  the  experimental  side,  LEED 
I-V  beam  profiles  have  been  measured  for  clean  surfaces  and  surfaces  with 
ordered  overlayers,  in  particular,  the  IrMll),  Ir(110),  Rh(lll)  and  Ni(lll) 
surfaces. 


The  formulation  of  a direct  method  of  analyzing  LEED  I-V  beam  profiles 
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to  arrive  at  surface  relaxation  via  a Fourier  transform  method  involved  the 
following  features:  (1)  The  inner  potential,  V , which  is  crucial  for  con- 
verting energy  to  momentum  transfer,  is  treated  as  a variable  to  be  deter- 
mined rather  than  as  a known  input  parameter.  (2)  Rather  than  deconvolute 
the  Patterson  function  (which  leads  to  instabilities),  a calculated  convo- 
lution product  is  compared  with  the  observed  Patterson  function.  (3)  The 
convolution  product  is  the  convolution  of  a "window"  function  and  a series 
of  delta  functions  representing  the  layer  positions.  The  layer  positions 
are  chosen  arbitrarily,  and  then  the  amplitudes  of  the  delta  functions  are 
obtained  by  performing  a least  squares  fit  of  the  convolution  product  to 
the  observed  Patterson  function.  (4)  The  mean  square  deviation  of  this 
fitted  convolution  product  with  the  observed  Patterson  function  is  calculated. 
(5)  Finally,  this  deviation  (called  the  residual)  is  plotted  as  a function 
of  both  Vq  and  the  degree  of  relaxation  (t)  of  the  surface  layer.  The  value 
of  t for  which  the  minimum  occurs  in  the  residual  gives  the  surface  structure 
in  the  I-V  profiles  in  the  sense  that  as  the  multiple  scattering  structure 
increases,  the  ability  to  determine  the  structure  diminishes  gradually,  not 
suddenly. 

This  method  has  been  illustrated  fully  using  "model"  kinematic  I-V  data 
from  a surface  which  is  both  relaxed  and  unrelaxed  {1) . It  was  found  that 
even  ideal  "data"  cannot  be  analyzed  to  give  an  "exact"  solution,  i.e.,  a 
perfect  fit  between  the  Patterson  function  of  the  kinematic  data  and  the 
convolution  product  of  the  window  transform  function  and  the  series  of  delta 
functions  representing  the  crystal  layer  spacing.  Nevertheless,  it  was 
shown  that  the  method  is  ctoable  of  determining  relaxations  even  from  data 
which  are  highly  dynamic.  This  was  illustrated  by  analyzing  the  theoretical 
I-V  profiles  of  the  (00),  (11),  (02)  and  (20)  beams  from  the  clean  W(110) 


surface  calculated  by  M.  A.  Van  Hove  and  S.  Y.  Tong  using  a multiple 
scattering  program  U,  2,  3).  The  correct  structure  was  determined  for 
surface  relaxations  of  + 9%,  OX  and  - 9X.  These  results  were  important  in 
that  they  showed  how  the  data  range  in  LEEO  affects  one's  ability  to  deter- 
mine surface  relaxation,  that  the  uncertainties  in  the  theory  are  not 
sufficient  to  destroy  the  ability  to  determine  the  structure,  and  that  the 
method  is  capable  of  resolving  relaxations  on  the  order  of  ± 5X  which  is 
comparable  to  the  uncertainty  in  multiple  scattering  analyses. 

Next,  the  transform-convolution  method  was  applied  to  experimental 

data  from  the  (111)  surface  of  Ir  measured  in  our  laboratory.  Multiple 

scattering  calculations  were  performed  also  to  determine  the  relaxation  of 

this  surface  for  the  first  time  (£,  5_).  Seven  specular  I-V  beam  profiles 

were  measured  from  15  to  975  eV  at  incident  angles  from  7°  to  62.5°  relative 

to  the  surface  normal.  The  transform-convolution  method  of  analysis  showed 
that  the  outermost  layer  of  Ir(lll)  is  either  unrelaxed  or  contracts  by  at  most 
4X  of  the  bulk  spacing.  In  agreement  with  this,  the  dynamic  calculations 
indicated  a contraction  of  2.5  ± 5%.  The  dynamic  calculations  showed  also  that 

the  registry  of  the  first  layer  of  this  crystal  surface  is  not  shifted,  main- 
taining the  fee  structure. 

Published  experimental  data  for  the  (110)  surface  of  the  fee  metals  Ni, 

A1  and  Ag  were  analyzed  using  the  transform-convolution  method  since  earlier 
dynamic  calculations  had  Indicated  a contraction  of  each  of  these  surfaces. 
The  results  of  the  transform-convolution  analysis  indicated  that  the  first 
layer  spacing  of  Ni (110) , Al(ilO)  and  Ag(110)  Is  contracted  by  5%,  4%  and  7% 
of  the  bulk  spacing,  respectively  (£,  6).  These  results  are  In  good  agree- 
ment with  the  results  which  had  been  obtained  previously  using  dynamic  calcu- 
lations. 


In  a similar  study,  both  specular  and  nonspecular  LEED  I-V  beam  profiles 
were  measured  for  the  Rh(lll)  surface  in  our  laboratory,  and  the  surface 
relaxation  was  determined  using  both  the  transform-convolution  method  and 
dynamic  calculations  (1).  Within  experimental  error,  neither  expansion  nor 
contraction  of  the  topmost  layer  was  detected.  The  results  of  the  transform- 
convolution  analysis  of  specular  beams  at  two  angles  of  incidence  and  of  a 
nonspecular  beam  at  normal  incidence  suggested  an  expansion  of  the  topmost 
layer  of  3 ± 5%  of  the  bulk  layer  spacing.  In  agreement  with  this,  com- 
parisons between  the  results  of  the  dynamical  calculation  and  the  experi- 
mental data  for  five  nonspecular  beams  at  normal  incidence  suggested  that 
the  surface  layer  relaxes  by  0 ± 5%.  In  addition,  the  dynamic  calculations 
indicated  that  the  topmost  layer  maintains  an  fee  structure. 

Two  other  clean  surface  structures  have  been  analyzed  using  multiple 
scattering  calculations,  the  Co(001)  surface  and  the  unreconstructed  Ir(110) 
surface.  The  surface  structure  of  Co(001)  was  examined  in  order  to  determine 
the  existence  of  residual  surface  phases  or  "frozen  in"  surface  embryos 
responsible  for  the  nucleation  of  the  Martensitic  transformation  in  Co  (8). 

It  was  found  that  no  such  residual  surface  embryos  exist,  and  that  the 
structural  phase  of  the  surface  is  equivalent  to  the  respective  allotropic 
phase  of  the  bulk.  Moreover,  the  transformation  at  the  surface  was  observed 
to  be  strongly  first-order  with  significant  temperature  hysteresis. 

An  Ir(110)  - (lxl)  surface  structure  was  prepared  by  adsorbing  a quarter- 
monolayer  of  oxygen  at  575  K on  a clean,  reconstructed  (1x2) surface  (9).  A 


The  experiments  were  performed  by  A.  Ignatiev  et  al.  at  the  University  of 
Houston,  and  the  calculations  were  performed  by  M.  A.  Van  Hove  at  Caltech. 
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comparison  between  measured  and  dynamically  calculated  LEED  I-V  beam  profiles 
Indicated  that  the  oxygen  is  distributed  randomly  over  the  crystal  surface, 
and  the  (lxl)  structure  is  the  same  as  a clean  unreconstructed  (lxl)  structure 

O 

with  a topmost  interlayer  Ir  spacing  of  1.26  ± 0.05  A.  This  is  equivalent  to 

O 

a contraction  of  approximately  7.5%  of  the  bulk  interlayer  spacing  of  1.36  A. 

Dynamic  calculations  were  performed  successfully  in  order  to  determine 
the  locations  of  hydrogen  atoms  on  a Ni (111)  surface  marking  the  first  deter- 
mination of  the  structure  of  a hydrogen  overlayer  on  any  surface  (10,  JU,  12) . 
Hydrogen  adsorption  below  room  temperature  on  Ni(lll)  yielded  the  appearance 
of  a (2x2)  overlayer  structure  of  which  the  "extra"  beams  attained  a maximum 
intensity  after  an  exposure  of  approximately  2.5  x 10'®  torr-sec.  A further 
Increase  in  the  exposure  led  to  a gradual  weakening  and  finally  complete 
disappearance  of  the  superstructure.  I-V  beam  profiles  were  measured  at 
maxiimm  intensity  at  a surface  temperature  of  110  K with  normal  incidence  of 
the  primary  beam.  These  beam  profiles  were  analyzed  by  dynamic  calculations 
using  reverse  scattering  perturbation  theory  in  conjunction  with  renormalized 
forward  scattering  and  layer  doublin  The  Ni  substrate  was  found  to  be 
unaffected  by  the  presence  of  the  hydrogen  atoms  which  are  arranged  in  a 
graphitic-like  overlayer  structure  with  a (2x2)  unit  cell,  i.e.,  both  types 
of  threefold  hollow  sites  on  this  fcc(lll)  surface  are  occupied.  No  detect- 
able buckling  of  the  hydrogen  overlayer  was  observed  to  occur.  The  overlayer- 
substrate  spacing  was  found  to  be  1.15  ± 0.1  X,  implying  a Ni-H  bond  length 
of  1.84  t 0.06  A. 

Finally,  dynamic  calculations  have  been  carried  out  to  determine  the 
location  of  the  oxygen  atoms  In  a c(2x2)  overlayer  structure  on  the  unre- 
constructed (110)  surface  of  Ir  (1^).  The  ordered  c(2x2)  oxygen  structure 
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was  prepared  by  chemisorbing  a ha If -monolayer  of  oxygen  at  room  temperature 
on  an  unreconstructed  (lxl)  Ir(110)  surface  stabilized  by  a quarter-mono’ ayer 
of  randomly  adsorbed  oxygen.  Results  of  the  structural  analysis  showed  that 
the  ordered  oxygen  atoms  are  residing  on  the  short-bridged  sites  (along  the 
rows)  on  the  (110)  surface.  The  ordered  overlayer-substrate  interlayer 

O 

spacing  was  found  to  be  1.37  ± 0.05  A,  corresponding  to  an  Ir-0  bond  length 

O 

of  1.93  ± 0.07  A.  The  topmost  substrate  interlayer  spacing  was  found  to  be 

o o 

1.33 ±0,07  A,  rather  than  1.26  ± 0.07  A which  is  the  topmost  interlayer 
spacing  of  the  unreconstructed  (lxl)  Ir(110)  surface. 

C.  Other  LEED  Results 

In  other  LEED  work,  not  related  to  determining  surface  structures,  th-ee 
contributions  have  been  made.  First,  a photographic  technique  has  been 
developed  which  allows  an  easy  and  accurate  determination  of  both  the  polar 
and  azimuthal  angles  of  incidence  in  a LEED  experiment  (1£).  This  method 
involves  ana’yzing  the  positions  of  the  diffraction  spots  on  a photograph  of 
the  LEED  pattern.  The  position  of  a single  spot  defines  uniquely  the  polar 
angle  (with  respect  to  the  surface  normal)  and  the  azimuthal  angle  (with 
respect  to  an  arbitrarily  defined  axis  in  the  plane  of  the  crystal  surface). 

A computer  program  has  been  developed  which  uses  as  input  data  the  beam 
energy,  the  beam  label  [the  value  of  (hk)],  the  spot  position  on  the  photo- 
graph, the  geometry  of  the  apparatus,  and  the  location  of  the  camera.  The 
program  corrects  for  distortions  caused  by  the  camera  and  gives  both  angles 
determined  by  each  spot  as  well  as  the  uncertainties  due  to  the  inability  to 
determine  the  spot  position  exactly.  For  a sinq'e  setting  of  the  crystal 


manipulator,  several  photographs  can  be  obtained  at  different  values  of  the 
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Incident  energy;  and  from  a single  photograph,  typically  on  the  order  of  ten 
spots  can  be  analyzed.  Using  this  large  number  of  angle  determinations, 
absolute  accuracies  better  than  0.1°  for  both  polar  and  azimuthal  incidence 
angles  are  obtained  routinely. 

Second,  the  Debye  temperature  of  the  Rh(lll)  surface  as  well  as  that  of 
bulk  Rh  has  been  measured  by  LEED  (7).  The  determination  of  the  Debye  temper- 
ature (or,  equivalently,  the  mean  square  displacement  of  the  atoms  due  to 
thermal  vibrations)  is  made  by  measuring  the  slope  of  the  logarithm  of  the 
intensity  of  the  specularly  diffracted  beam  as  a function  of  temperature. 

In  the  Debye  temperature  measurement  of  Rh(lll),  specular  LEED  beam  intensities 
at  several  incidence  angles  were  monitored  as  a function  of  temperature  over 
a range  of  electron  energies  from  approximately  30  eV  to  1000  eV.  It  was 
found  that  the  bulk  Debye  temperature  is  380  ± 23  K (determined  in  the  limit 
of  high  incident  beam  energy),  and  the  normal  component  of  the  Debye  temper- 
ature at  the  lowest  electron  energy  used  is  197  ± 12  K.  The  latter  may  be 
associated  with  the  normal  component  of  the  surface  Debye  temperature. 

Third,  a Monte  Carlo  simulation  has  been  carried  out  to  describe  two- 
dimensional  order-disorder  phenomena  measured  by  LEED  (15.,  16).  The  model 
contains  (attractive)  first,  (repulsive)  second,  and  (attractive)  third 
neighbor  pairwise  interactions.  The  special  case  of  oxygen  chemisorption  on 
a W(UO)  surface  iising  the  experimental  data  of  T.-M.  Lu,  G.-C.  Wang  and  M.  G. 
Lagally),  on  which  an  ordered  p(2xl)  overlayer  is  formed  atlowsurface  temperatures, 
was  considered  explicitly.  From  the  measured  order-disorder  transition 
temperatures  at  both  quarter-  and  ha If -mono layer  surface  coverages,  (non- 
unique) values  of  the  three  pairwise  interaction  energies  were  determined. 

These  pairwise  interaction  energies  were  used  to  determine  the  variation  in 
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the  total  interaction  energy,  the  heat  capacity  and  the  entropy  with  surface 
temperature. 

D.  Desorption  Phenomena 

Desorption  from  solid  surfaces  has  been  investigated  from  two  different 
points  of  view,  both  thermal  desorption  and  electron  beam  stimulated  desorp- 
tion. First,  a theoretical  analysis  of  thermal  desorption  mass  spectrometry 
(TDS)  was  performed  from  the  point  of  view  of  appropriate  dimensionless 
quantities  (17).  Analytical  expressions  were  derived  which  allow  a deter- 
mination of  the  activation  energy  of  desorption  and  the  pre-exponential 
factor  of  the  desorption  rate  coefficient  using  parameters  obtained  easily 
from  thermal  desorption  mass  spectra.  In  particular,  it  was  demonstrated 
that  the  spectral  peak  widths  and  the  temperature  at  which  the  maximum  rate 
of  desorption  occurs  may  be  used  to  describe  both  first  and  second  order 
desorption  kinetics.  The  use  of  the  method  was  illustrated  explicitly  by 
applying  it  to  the  analysis  of  several  important  classes  of  desorption 
reactions.  Then,  the  distortions  which  are  associated  with  TDS  as  a result 
of  the  combined  effect  of  surface  heating  rate  and  system  pumping  speed  were 
quantified  theoretically  (18).  It  was  demonstrated  that  considerable  error 
In  the  determination  of  the  activation  energy  and  the  pre-exponential  factor 
of  desorption  can  be  introduced  by  the  combined  effects  of  high  heating  rate 
and  low  pumping  speed.  It  was  found  that  for  an  accurate  determination  of 
the  activation  energy  and  the  pre-exponential  factor  for  first-order  desorp- 
tion, the  reciprocal  of  the  product  of  the  heating  rate  and  the  pumping  time 
constant  must  be  rather  large,  e.g.,  greater  than  approximately  0.5. 

The  other  aspect  of  desorption  which  was  considered  is  the  case  of  ESO. 
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Since  the  electron  beam  used  as  a probe  in  electron  scattering  experiments 

such  as  LEEO  and  Auger  electron  spectroscopy  (AES)  may  perturb  a chemisorbed 

overlayer,  it  is  imperative  to  investigate  this  effect.  For  this  reason, 

electron  beam  induced  perturbations  of  CO  chemisorbed  on  Ir(lll)  were 

Measured  experimentally  using  LEED  and  AES  (19).  The  total  interaction 

cross  section  for  electron  stimulated  desorption  and  dissociation  was  found 

-17  2 

to  be  approximately  0.8  to  1.7  x 10  cm  when  defined  with  respect  to  the 
primary  flux  of  a 2.5  keV  electron  beam.  Electron  stimulated  dissociation 
was  found  to  occur  at  1-21  of  the  rate  of  electron  stimulated  desorp- 
tion. 

E.  Model  Calculations:  Inelastic  Electron  Tunneling  Spectroscopy, 
Surface  Reconstruction,  Changes  in  Electronic  Densities  of  States 

A theory  was  developed  which  describes  the  dependence  of  the  spectral 
intensity  upon  adsorbate  concentration  in  inelastic  electron  tunneling 
spectroscopy  (20-22).  In  contrast  to  previous  theories  which  predicted  a 
linear  dependence,  it  was  shown  that  the  intensity  of  molecular  vibrational 
loss  peaks  varies  approximately  as  n*'\  where  n fs  the  concentration  of 
adsorbate  on  the  Insulator  surface  in  a metal-i'isulator-metal  tunneling 
junction.  The  actual  dependence  does  not  follow  a simple  power  law,  but 
over  the  range  of  from  51  to  901  of  saturation  coverage,  the  dependence  can 
be  approximated  by  n1'^.  This  result  Is  in  agreement  with  experimental 
Measurements  of  J.  Langan  and  P.  Hansma  who  examined  benzoic  acid  chemi- 


sorption on  an  aluminum  oxide  surface.  This  model  makes  It  possible  to  use 
intensity  measurements  to  determine  relative  coverages  on  the  surface,  and 
consequently,  evaluate  the  kinetics  of  adsorption. 
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A Green's  function  perturbation  technique  was  developed  which  is 
appropriate  for  obtaining  the  change  in  the  total  electronic  energy  when 
the  surface  atoms  of  a crystal  reconstruct  in  a periodic  manner  (23).  The 
phase  shift  technique  was  used  to  determine  the  change  in  the  electronic 
density  of  states.  The  method  is  quite  general  and  can  be  used  to  study 
arbitrary  surfaces  which  are  reconstructed.  As  a specific  example,  the 
change  in  the  total  electronic  energy  which  occurs  upon  reconstruction  of 
the  (001)  surface  of  a model  two-band  crystal  with  the  CsCl  structure 
was  calculated  (24).  The  energy  changes  due  to  a (2x1)  and  a c(2x2) 
reconstruction  were  considered.  The  occurrence  of  reconstruction  on  a surface 
was  shown  to  be  dependent  upon  the  relative  magnitudes  of  the  electronic 
energy  and  the  elastic  strain  energy. 

Finally,  a model  calculation  was  carried  out  for  the  orbital  density 
of  states,  as  well  as  the  change  in  the  density  of  states  due  to  the  chemi- 
sorption, of  a two-level  adsorbate  bonded  to  the  s-band  of  an  fcc(lOO)  metal 
surface  (25).  The  calculation  was  carried  out  for  an  adsorption  geometry  with  the  adatom 
situated  over  the  four-fold  hole  site,  with  a n-bonding  interaction  with  the  diagonal 
substrate  atoms.  A Green's  function  formalism  was  used  within  the  LCAO-tight 
binding  approximation.  It  was  found  that  orbital  resonances  contain  not  only 
contributions  from  the  adorbitals  but  also  from  the  substrate  group  orbital 
which  participates  in  the  bonding.  The  admixture  of  each  orbital  in  the 
resonances  can  be  understood  qualitatively  in  terms  of  both  direct  and 
indirect  interactions  which  depend  on  the  parameters  of  the  model,  namely, 
the  unperturbed  adorbital  energies,  the  adsorbate-substrate  coupling  strengths, 
and  the  intra-adsorbate  coupling  strength. 
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F.  Chemisorption  and  Surface  Reactions 

In  addition  to  surface  structural  determinations,  recent  efforts  have  focused 
strongly  on  investigations  of  chemisorption  and  surface  reaction  kinetics.  This 
work  will  continue  with  the  support  of  a renewal  proposal  now  under  review  by  the 
ARO.  Work  already  completed  in  this  area  includes  the  refinement  of  a method  of 
measuring  work  function  changes  continuously  and  a study  of  the  chemisorption  of 
CO  on  the  (110)  surface  of  Ir. 

First,  a method  was  developed  for  measuring  work  function  changes  con- 
tinuously (26).  It  is  an  ac  variation  of  the  diode  technique  and  can  be  imple- 
mented easily  in  a standard  retarding  potential  LEED/Auger  system  with 
derivative  detection,  ileasurements  for  oxygen  adsorption  on  the  (110)  surface 
of  Ir  were  made  to  illustrate  that  this  method  can  resolve  work  function 
changes  within  5 meV. 

Also,  the  chemisorption  of  CO  on  well-characterized  clean  and  oxidized 

Ir(110)  surfaces  was  studied  in  an  ultrahigh  vacuum  environment  (27_) . On  both 

surfaces,  it  was  found  that  adsorption  occurs  via  a mobile  precursor  state. 

At  surface  temperatures  of  both  90  and  300  K,  the  adsorption  kinetics  were 

shown  to  be  independent  of  whether  the  surface  is  clean  or  oxidized.  Both 

LEED  and  thermal  desorption  data  indicated  that  saturation  coverage  is 
15  2 

approximately  10  molecules/cm  , or  one  monolayer  under  all  conditions.  No 
evidence  for  physical  adsorption  of  CO  was  found,  even  at  90  K,  at  which 
temperature  the  probability  of  adsorption  is  nearly  unity  until  saturation 
coverage  is  reached.  At  300  K,  the  probability  of  adsorption  was  found  to 
be  unity  at  low  surface  coverages  but  decreased  dramatically  at  0.8  of 
saturation  coverage.  This  decrease  was  attributed  to  interadsorbate  repulsion. 

A kinetic  model  was  formulated  which  suggests  that  a repulsive  interaction  of 
2 kcal/mole  exists  between  CO  molecules  at  that  fractional  surface  coverage. 
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T'*  this  point,  that  work  which  has  been  completed  has  been  summarized. 
In  addition,  a significant  amount  of  experimental  information  has  been 
■easured  for  the  adsorption  of  H^,  CO  and  the  nonreactive  coadsorption 
of  and  0^,  and  and  CO;  and  the  reactive  coadsorption  of  and  O2  on 
the  (111)  surface  of  Rh.  Results  to  date  will  be  presented  below,  although 
this  work  will  be  completed,  and  manuscripts  will  be  prepared  with  the 
support  of  the  renewal  grant. 


1.  The  Adsorption  of  Oxygen  on  the  (111)  Surface  of  Rhodium,  P.  A.  Thiel, 

J.  T.  Yates,  Jr.  and  W.  H.  Weinberg  (in  preparation). 

The  adsorption  of  oxygen  on  Rh(lll)  at  100  K has  been  studied  by  TDS, 

AES  and  LEED.  Oxygen  adsorbs  into  a disordered  state  at  100  K and  orders 
irreversibly  into  an  apparent  (2x2)  surface  structure  upon  heating  to 
temperatures  above  150  K.  On  the  (111)  surface,  this  could  either  be  a true 
(2x2)  overlayer  or  three  independent  domains  of  (2x1)  structures  rotated 
120°  with  respect  to  one  another.  LEED  I-V  data  have  been  collected,  and  a 
dynamic  LEED  calculation  will  be  performed  to  clarify  the  structure  of  the 
"(2x2)"  oxygen  overlayer.  The  kinetics  of  ordering  of  the  oxygen  have 
been  measured  at  surface  temperatures  between  150  and  280  K by  monitoring 
the  intensity  of  the  (l,*s)  LEED  beam  as  a function  of  time  at  various  surface 
temperatures.  The  kinetic  data  are  consistent  only  with  a model  in  which 
the  rate  of  ordering  of  oxygen  adatoms  is  proportional  to  the  square  of  the 
concentration  of  disordered  oxygen.  The  activation  energy  for  ordering  was 
found  to  be  13.5  ± 0.5  kcal/mole.  The  fact  that  the  ordering  kinetics  are 
proportional  to  the  square  of  the  concentration  of  disordered  oxygen  may 


■p 
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well  be  due  to  the  nature  of  the  interactions  among  the  oxygen  adatoms  on 
the  surface,  e.g.,  t le  additionof  a single  adatom  to  the  ordered  structure 
nay  be  energetically  unfavorable  compared  to  the  concerted  addition  of  two 

■ 

adatoms  to  the  ordered  structure.  This  point  of  view  will  be  formulated 
quantitatively  in  an  effort  to  explain  the  unexpected  second  order  kinetics 
which  were  observed.  At  surface  temperatures  above  280  K,  the  oxygen  under- 
goes a two-step  irreversible  disordering  and  dissolution  process.  Formation 
of  the  high  temperature  disordered  state  is  retarded  at  high  coverages  of 
oxygen  on  the  surface. 

The  TDS  and  AES  data  show  also  that  oxygen  adlayers  are  depleted  by 
dissolution  into  the  bulk  of  the  Rh  at  temperatures  above  400  K.  The  thermal 
desorption  of  oxygen  shows  a "threshold  effect",  with  no  desorption  of 
oxygen  occurring  below  an  exposure  of  approximately  1.3  x 10~6  torr-sec,  even 
though  a significant  quantity  of  oxygen  chemisorbs  at  lower  exposures  as 
judged  by  AES.  Analysis  of  the  thermal  desorption  data  at  rather  low  surface 
coverages  (extrapolated  to  zero  coverage)  yields  values  of  the  pre-exponential 
factor  of  the  (second  order)  desorption  rate  coefficient  and  the  activation 
energy  of  desorption  of  2.5  x 10  cm  /sec  and  56  ± 2 kcal/mole,  respectively. 
At  higher  coverages,  the  desorption  data  are  complicated  by  contributions 
from  multiple  states.  At  100  K,  the  probability  of  adsorption  of  oxygei  is 

' 

approximately  0.2,  a value  which  is  maintained  to  rather  high  surface 
Coverages  due  to  the  precursor  nature  of  the  adsorption  kinetics. 

2.  The  Adsorption  of  Hydrogen  on  the  (111)  Surface  of  Rhodium,  J.  T.  Yates, 
Jr.,  P.  A.  Thiel  and  W.  H.  Weinberg  (in  preparation). 


The  adsorption  and  desorption  kinetics  of  hydrogen  on  Rh(lll)  have  been 


studied  using  TDS.  The  adsorption  of  hydrogen  at  temperatures  below  175  K 
occurs  with  a probability  of  approximately  0.2.  At  low  surface  coverages, 
a desorption  state,  exhibiting  second-order  desorption  kinetics,  is  observed 
near  375  K.  As  the  coverage  of  hydrogen  increases,  the  desorption  peak 
shifts  to  lower  temperatures  and  broadens  due  to  interact-'onal  effects  in 
the  hydrogen  overlayer.  No  ordering  of  the  hydrogen  on  the  surface  is 
observed  by  LEED  at  any  coverage.  Following  adsorption  of  a mixture  of 
hydrogen  and  deuterium,  isotopic  mixing  is  observed  to  produce  ^ and 
HD.  The  desorption  of  pure  deuterium  was  observed  to  occur  with  identical 
kinetic  parameters  to  pure  hydrogen.  A kinetic  model  will  be  formulated 
describing  both  the  adsorption  and  the  desorption  kinetics.  In  the  case  of 
the  latter,  the  variation  with  surface  coverage  of  both  the  pre-exponential 
factor  of  the  desorption  rate  coefficient  as  well  as  the  activation  energy 
of  desorption  w.  1 be  evaluated. 

3.  The  Adsorption  of  Carbon  Monoxide  on  the  (111)  Surface  of  Rhodium, 

E.  D.  Williams,  P.  A.  Thiel,  J.  T.  Yates,  Jr.  and  W.  H.  Weinberg  (in 

preparation). 

The  adsorption  of  CO  on  Rh(lll)  at  90  K has  been  studied  by  TDS  and  LEED. 
The  LEED  data  suggest  that  CO  adsorbs  in  a number  of  ordered  structures  which 
are  produced  at  different  coverages  on  the  surface.  The  progression  of  CO 
LEED  patterns  observed,  from  low  to  high  coverages,  is  the  following:  (1)  A 
(/T  x /T  )R  30°  structure  at  fractional  CO  coverages  (0^)  below  a third  of 
a monolayer,  (2)  Various  "split"  (2x2)  structures  at  6^  > 1/3,  and  (3)  An 
ordered  (2x2)  structure  at  saturation  coverage.  This  progression  indicates 
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a continuous  compression  of  the  CO  overlayer  with  increasing  CO  coverage,  a 
phenomenon  which  has  been  observed  previously  in  our  laboratory  in  the  case 
of  CO  chemisorption  on  Ir(lll)  and  Ru(001).  The  relationship  between 
intensity  and  exposure  of  both  the  (/ T x /3)R  30°  and  the  (2x2)  LEED 
structures  have  been  measured  and  will  be  related  to  the  coverage-exposure 
functional  measured  with  TDS. 

Thermal  desorption  spectra  of  CO  from'Rh(lll)  yield  a first-order 
desorption  feature  which  shifts  to  lower  temperature  and  broadens  as  coverage 
increases,  a result  of  interadsorbate  interactions.  The  kinetic  parameters 
of  the  desorption  rate  coefficient  will  be  evaluated  at  low  coverages,  and, 
if  possible,  as  a function  of  surface  coverage.  An  upper  limit  for  the  ESD 
cross  section  of  CO  on  Rh(lll)  will  be  evaluated  from  the  measured  LEED  data. 
This  Is  of  importance  to  latter  work  in  which  the  structure  of  CO  in  the 
(/J  x /3)R  30°  overlayer  at  0gg  = 1/3  will  be  determined  by  LEED. 

4.  The  Coadsorption  of  Hydrogen  and  Oxygen  on  the  (111)  Surface  of  Rhodium 

at  90  K,  P.  A.  Thiel,  J.  T.  Yates,  Jr.  and  W.  H.  Weinberg  (in  preparation). 

The  coadsorption  of  oxygen  and  hydrogen  on  Rh(lll)  has  been  studied  by 
TDS,  AES  and  LEED.  An  important  finding  in  this  work  is  that  the  coadsorption 
phenomena  which  were  observed  depend  strongly  upon  which  of  the  components  of 
the  coadsorbed  overlayer  is  chemisorbed  initially,  i.e.,  there  is  a sensitive 
dependence  upon  the  order  of  adsorption. 

First,  the  coverage  of  hydrogen  on  Rh(lll)  has  been  studied  both  as  a 
function  of  oxygen  exposure  at  90  K prior  to  hydrogen  exposure  at  90  K,  and 
as  a function  of  the  temperature  to  which  a saturated  oxygen  overlayer  Is 
annealed  prior  to  hydrogen  exposure  at  90  K.  Exposing  a clean  Rh(lll) 
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surface  to  oxygen  at  90  K prior  to  exposing  it  to  hydrogen  at  the  same 
temperature  partially  blocks  hydrogen  adsorption  and  causes  the  hydrogen  to 
desorb  in  two  new  states  at  approximately  200  K and  400  K.  However, 
annealing  the  oxygen  overlayer  allows  the  adsorption  of  up  to  50%  of  the 
saturation  amount  of  hydrogen  on  clean  Rh(lll)  as  judged  by  IDS.  This  proves 
that  hydrogen  adsorption  and  desorption  are  extremely  sensitive  to  the 
disorder-order-disorder  transitions  of  the  oxygen  overlayer  discussed  earlier 
(Sect.  F-l). 

Second,  the  adsorption  of  oxygen  on  a Rh(lll)  surface  pre-exposed  to 
hydrogen  was  investigated.  When  the  Rh(lll)  surface  is  saturated  with 
hydrogen  at  90  K and  then  given  a saturation  exposure  of  oxygen  also  at  90  K, 
the  hydrogen  desorbs  within  10  K of  the  peak  temperature  of  desorption  on  the 
clean  surface,  with  almost  no  attenuation  of  the  intensity  of  the  desorption 
peak.  This  result  is  surprising  in  view  of  the  fact  that  AES  indicates  that 
the  amount  of  oxygen  adsorbed  on  the  Rh(lll)  surface  saturated  with  hydrogen 
Is  attenuated  by  only  20%  relative  to  the  clean  Rh(lll)  surface.  However, 
no  oxygen  thermal  desorption  occurs  from  this  coadsorbed  overlayer,  in  marked 
contrast  to  the  behavior  of  oxygen  on  clean  Rh(lll)  as  discussed  previously 
in  Section  F-l.  The  thermal  desorption  of  oxygen  from  the  surface  pre-exposed 
to  hydrogen  has  been  monitored  by  TDS  as  a function  of  hydrogen  surface 
coverage.  The  pre-adsorbed  hydrogen  blocks  oxygen  desorption  completely  for 
eH  i 0.6,  where  0^  is  the  fractional  coverage  of  hydrogen  relative  to  satur- 
ation coverage.  Moreover,  LEED  measurements  indicate  that  pre-adsorbed 
hydrogen  prevents  oxygen  from  ordering  on  the  surface  at  temperatures  up  to 
the  onset  of  oxygen  desorption,  800  K,  a temperature  at  which  no  hydrogen 
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remains  or  the  surface.  The  presence  of  disordered  oxygen  on  the  surface 
was  verified  by  AES,  and  TDS  was  used  to  show  that  no  detectable  water  was 
produced  under  these  experimental  conditions.  The  production  of  water  from 
hydrogen  and  oxygen  on  Rh(lll)  is  discussed  below  in  Section  F-6.  Conse- 
quently, it  may  be  concluded  that  the  presence  of  hydrogen  on  Rh(lll) 
influences  profoundly  the  subsequent  adsorption  of  oxygen,  even  though  oxygen 
desorption  occurs  at  temperatures  far  higher  than  those  at  which  hydrogen 
desorption  occurs.  The  absence  of  any  oxygen  thermal  desorption  from  a 
Rh(lll)  surface  pre-exposed  to  large  amounts  of  hydrogen  may  be  related  to 
the  lack  of  formation  of  the  ordered  (2x2)  superstructure  which  occurs  on 
the  clean  Rh(lll)  surface.  Additional  work  relating  to  this  nonreactive 
coadsorption  system  will  be  carried  out  during  the  initial  stages  of  support 
of  the  renewal  grant. 

5.  The  Cuadsorption  of  Carbon  Monoxide  and  Hydrogen  on  the  (111)  Surface 

of  Rhodium,  P.  A.  Thiel,  J.  T.  Yates,  Jr.,  E.  D.  Williams  and  W.  H. 

Weinberg  (in  preparation). 

• • 

Both  TDS  and  LEED  have  been  used  to  study  the  nonreactive  interaction 
of  hydrogen  and  CO  on  Rh(lll).  When  a Rh(lll)  surface  saturated  with 
hydrogen  is  exposed  to  CO  at  125  K,  no  displacement  of  hydrogen  is  observed. 
However,  the  peak  in  the  hydrogen  thermal  desorption  spectrum  shifts  mono- 
tonlcally  to  lower  temperatures  as  a function  of  CO  coverage,  indicating  a 
decrease  In  the  activation  energy  for  hydrogen  desorption.  In  contrast  to 
this  behavior,  preadsorption  of  CO  at  temperatures  below  273  K,  followed  by 
a saturation  exposure  of  hydrogen  at  125  K,  leads  to  an  attenuation  of  the 
Intensity  of  the  hydrogen  thermal  desorption  spectrum.  However,  this 
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attenuation  is  not  accompanied  by  a shift  in  the  peak  temperature  indicating 
that  the  activation  energy  for  desorption  is  not  changed  relative  to  that 
which  obtains  for  the  clean  Rh(lll)  surface. 

As  discussed  in  Section  F-3,  exposure  of  the  clean  Rh(lll)  surface  to 
CO  at  100  K leads  to  a succession  of  ordered  structures  indicative  of  a 
continuous  compression  of  the  CO  overlayer  with  increasing  coverage.  The 
first  ordered  structure,  a (/ J x /3)R  30°  structure,  exhibits  a maximum 
LEED  intensity  after  an  exposure  of  approximately  1.5  x 10"®  torr-sec  of  CO 
onto  the  initially  clean  surface.  Preadsorption  of  hydrogen  causes  this 
intensity  maximum  to  decrease  and  to  shift  to  lower  CO  exposures.  Similarly, 
post-exposure  of  hydrogen  to  a surface  on  which  a well  ordered  (/3  x /J)R  30° 
CO  overlayer  exists  causes  the  Intensity  of  the  CO  superstructure  to  decrease. 

As  expected,  no  methane  desorption  was  observed  from  the  coadsorbed 
hydrogen  and  CO  overlayers  at  a sensitivity  level  of  < 0.5%  of  desorption  of 
a monolayer  of  pure  CO. 

6.  The  Catalytic  Reaction  between  Adsorbed  Oxygen  and  Hydrogen  on  the  (111) 
Surface  of  Rhodium,  J.  T.  Yates,  Jr.,  P.  A.  Thiel  and  W.  H.  Weinberg  (in 
preparation). 

Oxygen  was  found  to  chemisorb  into  a disordered  state  on  the  (111)  surface 
of  Rh  near  100  K.  As  discussed  in  Section  F-l,  an  activated  conversion  to  an 
ordered,  apparent  (2x2)  overlayer  structure  occurs  at  surface  temperatures 
between  approximately  150  and  280  K.  Above  280  K,  there  is  a second  activated 
process  with  an  activation  energy  of  8.2  ± 0.3  kcal/mole  which  leads  to  dis- 
ordered oxygen.  The  adsorption  of  oxygen  on  Rh(lll)  at  335  K occurs  via 
Langmuir  kinetics  with  a probability  of  adsorption  of  approximately  0.6. 


During  adsorption  at  this  temperature,  a mixture  of  ordered  and  disordered 

oxygen  is  produced  with  massive  conversion  to  the  ordered  layer  in  the  final 

stages  of  adsorption  at  exposures  of  oxygen  greater  than  5 x 10'6  torr-sec. 

At  fractional  surface  coverages  below  half  of  saturation,  the  intensity  of 

the  (1,5)  LEED  beam  increases  as  the  square  of  the  fractional  surface  coverage 

as  expected  for  the  formation  of  ordered  oxygen  islands. 

The  reaction  of  the  ordered  oxygen  overlayer  with  hydrogen  to  produce 
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water  occurs  rapidly  above  275  K with  hydrogen  pressures  in  the  10  to  10 

torr  region.  The  reaction  exhibits  first  order  kinetics  in  hydrogen  pressure. 

Implying  a trapping  mechanism  for  adsorbed  hydrogen  on  the  ordered  oxygen 

overlayer,  i.e.,  irreversible  adsorption  of  hydrogen.  The  activation  energy 

for  the  removal  of  ordered  oxygen  by  reaction  with  hydrogen  is  5.3  kcal/mole, 

°2 

and  the  reaction  cross  section  is  approximately  250  A (the  pre-exponential 
factor  of  the  reaction  rate  coefficient).  This  large  cross  section  for  the 
loss  of  ordered  oxygen  indicates  that  near  saturation  coverage  of  oxygen, 
massive  disordering  occurs  for  small  decreases  in  oxygen  coverage  during 
reaction.  This  observation  is  consistent  with  the  earlier  observation  of 
massive  oxygen  ordering  during  the  final  stages  of  oxygen  adsorption  at  335  K. 
A Langmuir-Hinshelwood  mechanism  for  the  reaction  is  suggested  by  the  obser- 
vation of  H2O,  HDO  and  D2O  reaction  products  when  mixtures  of  gaseous  hydrogen 
and  deuterium  were  used  to  react  with  the  chemisorbed  oxygen. 


1. 
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